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B-mode of polarized anisotropies of the cosmic microwave background is a unique and nearly 
direct probe of primordial inflation, which can constrain the amplitude of the primordial gravity 
waves. However, its detection and precise measurement is made difficult by a minute amplitude of 
the signal, which has to be discerned from many contributions of non-cosmological origin and reliable 
estimated in the presence of numerous sources of statistical uncertainties. Among these latter, the E- 
to-B leakage, arising as a result of partial sky coverage, has been found to play a key and potentially 
fundamental role in determining the possible statistical significance with which the primordial B- 
mode signal can be detected. In this work we employ the pure-pseudo formalism devised to minimise 
the effects of the leakage on the variance of power spectrum estimates and discuss the limits on the 
tensor-to-scalar ratio, r, that could be realistically set by current and forthcoming measurements 
of the B-mode angular power spectrum. We compare those with the results obtained using other 
approaches: naive mode-counting, minimum-variance quadratic estimators, and re-visit the question 
of optimizing the sky coverage of small-scale, suborbital experiments in order to maximize the 
statistical significance of the detection of r. We show that the optimized sky coverage is largely 
insensitive to the adopted approach at least for reasonably compact sky patches. We find, however, 
that the mode-counting overestimates the detection significance by a factor ~ 1.17 as compared to 
the lossless maximum variance approach and by a factor ~ 1.25 as compared to the lossy pure pseudo¬ 
spectrum estimator. In a second time, we consider more realistic experimental configurations. With 
a pure pseudospectrum reconstruction of B-modes and considering only statistical uncertainties, we 
find that a detection of r ~ 0.11, r ~ 0.0051 and r ~ 0.0026 at 99% of confidence level is within the 
reach of current sub-orbital experiments, future arrays of ground-based telescopes and a satellite 
mission, respectively. This means that an array of telescopes could be sufficient to discriminate 
between large- and small-field models of inflation, even if the E-to-B leakage is consistently included 
but accounted for in the analysis. However, a satellite mission will be required to distinguish between 
different small-field models depending on the number of e-folds. 
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I. INTRODUCTION 

Primordial gravity waves are expected to be produced 
during cosmic inflation in addition to scalar perturba¬ 
tions. If indeed present, they would leave a characteris¬ 
tic footprint on the polarized anisotropies of the cosmic 
microwave background (CMB), as they are considered 
to be essentially the sole source of the so-called primor¬ 
dial B-mode residing at the super-horizon scales at the 
time of the last scattering. A detection of the B-mode 
angular power spectrum at large angular scales would be 
then treated as a smoking gun of inflation, while a precise 
measurement of its amplitude would constrain the energy 
scale of inflation, or, geometrically speaking, the expan¬ 
sion rate of the Universe during inflation [TJ 2] ■ This 
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amplitude is expressed by the tensor-to-scalar ratio, r, 
defined as the relative power of primordial gravity waves 
with respect to that of the scalar perturbations at some 
pivot scale kg, chosen here to be equal to kg = 0.002 
Mpc~ 1 . Currently, the most stringent upper bound on r 
using temperature anisotropies has been derived by the 
Planck collaboration : r < 0.11 at 95% CL [3], while a 
recent joint analysis of the Planck and BICEP2 polarized 
data set an upper limit ro .05 < 0.12 at 95% CL [6j. The 
measurement of the tensor-to-scalar ratio r could allow 
to discriminate between different inflationary models. In 
particular, if this upper bound r ~ 0.1 is indeed real¬ 
ized in nature, this would imply a rather high amount of 
primordial gravity waves thus favoring large-fields infla¬ 
tionary models [?]. 

At smaller angular scales, the B-mode is dominated 
by the lensing induced signal. This signal is generated 
by the gravitational lensing of the CMB photons due to 
the large scale structure [8j. The lensing contribution 
is well-understood from a theoretical point of view and 
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can be uniquely predicted given the primary -E-modes 
anisotropies and the lensing deflection field [9j. Such 
predictions have been been recently confirmed by the 
SPTpol [ T01IIlj and POLARBEAR experiments dna, 
with also constraints on the CMB lensing E-mode power 
spectrum mm- The lensing E-mode signal does not 
depend on r. It therefore acts as a source of an additional 
’noise’ masking the primordial, r-dependent E-mode, 
and making its detection more difficult. Striving for 
a detection of r, one has to either try to remove this 
lensing signal [TB] or rely solely on the large angular 
scales. In this latter case two features of the primordial 
E-mode spectrum are of particular interest as they are 
anticipated to be particularly prominent. These are 
so-called reionization and recombination bumps peaking 
at 2 < i < 10 and at £ ~ 100, respectively. 

Measuring E-mode is made even more difficult by the 
fact that measurements as performed by the majority 
of current experiments, which scan the sky in order to 
produce its maps, are straightforwardly expressed only 
in terms of the Stokes parameters, Q and U. The E- 
and E-mode are mathematically related to the Stokes pa¬ 
rameters UzHIS] and can be therefore recovered from the 
observational data. This however is only simple, if full- 
sky data were available. In contrast, realistic CMB ex¬ 
periments provide maps of polarized anisotropies, which 
only cover a reduced fraction of the celestial sphere, rang¬ 
ing from ~ 1% for balloon-borne and ground-based ex¬ 
periments to ~ 70% for satellite missions. In the con¬ 
text of pseudospectrum estimation of the angular power 
spectra on an incomplete sky part of the E-mode signal 
is unavoidably mislabelled as E-modes and vice verse. 
Though such leakages can be corrected on average, the 
leaked signal inevitably contributes to the sampling vari¬ 
ance of the other reconstructed spectrum. This dramati¬ 
cally increases the uncertainties of the estimated E-mode 
spectrum since the cosmological E-mode is expected to 
be at least two orders of magnitude higher than the E- 
mode in terms of their power spectrum [20] . The na¬ 
ture of the leakages and approaches to their removal 
were investigated in Ref. m and a relevant pseudospec¬ 
trum estimator, referred to as the pure pseudospec¬ 
trum estimator, was proposed subsequently in Ref. |22|. 
This estimator has been thoroughly investigated and ex¬ 
tended to include an optimization of the sky apodiza- 
tion [52j, cross-spectrum approaches |2H, and, TB and 
EE cross-correlations [25j. Alternative constructions of 
pseudospectrum estimators correcting for E-to-E leak¬ 
ages have been also proposed [2B1 - E3] . Nevertheless, the 
pure pseudospectrum method has been found the most 
mature and efficient one, particularly due to its ability of 
optimizing the sky apodizations [30] , making it a method 
of choice for many practical applications. It is worth 
pointing out that the leakages are indeed ubiquities and 
correcting for them is as mandatory for small-scale ex¬ 
periments, covering ~ 1% of the sky, as for satellite-like 
missions, with access to as much as ~ 70% of the sky 3()| . 


Though the impact of the E-to-E leakage on the 
variance of the E-mode power spectrum is gener¬ 
ally acknowledged, it is rarely included in projecting 
performance of planned CMB experiments or instru¬ 
mental concepts from the point of view of their setting 
constraints on the tensor-to-scalar ratio, r. Instead, 
the major body of work (see EH EH for some recent 
examples) in this area is based on simplified mode¬ 
counting arguments (see, however e.g., [531ITTTT] for some 
exceptions). This stemmed mostly from the practical 
reasons, as the impact of the leakage is neither calculable 
analytically nor analysis method independent. 

The objective of this work is to fill this gap and present 
a more systematic study of the impact of the presence of 
the leakage on the performance forecasts of CMB B-mode 
experiments. The paper consists of two parts. In the 
first part, Sec. [TTJ we consider idealized observations of 
azimuthally symmetric sky areas with homogenous noise 
and study differences between performance forecasts de¬ 
rived applying three different approaches for different as¬ 
sumed sky area sizes. Subsequently, from these three 
different perspectives we revisit the issue of the optimal 
sky area, which would permit setting the most stringent 
constraints on the scalar-to-tensor ratio, r, given a fixed 
length and sensitivity of the experiment. 

In the second part, Sec. |III[ we complete those con¬ 
siderations by discussing more realistic sky areas defined 
for three types of experiments: small-scale observations 
covering ~ 1% of the sky, an array of ground-based tele¬ 
scope covering ~ 36% and a satellite-like mission capable 
of delivering up to 71% of the foreground clean sky. Our 
conclusions are drawn out in Sec. m where we also 
briefly sketch the implications for constraining inflation¬ 
ary models. 

Throughout this work we neglect complications such as 
polarized diffuse foregrounds, e.g., [331 Ml and account 
for resolved points sources only by appropriately tailoring 
the adopted mask. We also assume that no subtraction 
of the lensing E-mode has been attempted [IBl- 


ii. MEASURING THE TENSOR-TO-SCALAR 
RATIO FOR IDEALIZED SMALL-SCALE 
EXPERIMENT 

A. Experimental setup 

We consider first the case of small-scale experiments 
in an idealized way. The observed part of the celestial 
sphere is assumed to be azimuthally symmetric, given 
by a spherical cap. We however let vary the sky cover¬ 
age from 0.5% to 10%. The noise is an homogeneous, 
white noise, and its level is fixed at np(l%) = 5.75^iK- 
arcminute for f s k y = 1% (a typical level for ongoing 
small-scale experiments). For a fixed sensitivity and 
a fixed time of observation, the noise level (in /iK- 
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arcminute) scales as: 


nH/sky) = /^xn P ( 1%). (1) 

The instrumental noise reprojected on the sky thus varies 
from 4.1/iK-arcminute to 18/iK-arcminute for an ob¬ 
served fraction of the sky of 0.5% and 10%, respectively. 
Finally, the angular resolution is given by an azimuthally 
symmetric, gaussian beam with a width of 8 arminutes. 

We subsequently investigate the signal-to-noise ratio, 
(S/N) r = r/a y, as a function of the sky coverage. This 
will be done considering four values of the tensor-to- 
scalar ratio: r = 0.07, 0.1, 0.15 and 0.2 . We note 
that the last two values are disfavored by the current 
data, nevertheless we include them in our considerations 
as they are useful in demonstrating some of the effects 
we describe hereafter. 


B. Fisher matrix formalism 


Translating the uncertainties on the B-mode angular 
power spectrum reconstruction into error bars on the 
measured tensor-to-scalar ratio, ay, can be done using 
a Fisher matrix formalism. 

For the rather small observed fractions of the celes¬ 
tial sphere here-considered, the B-mode angular power 
spectrum is reconstructed within bandpowers, labelled b 
hereafter, with bandwidths Aj,. The binned power spec¬ 
trum is given by Cf = Ye Pm Cf , where the binning 
operator is defined as: 


Pm 


1 ( 1 + 1 ) 

27tA 6 

0 


if £ G b, 
if l^b. 


( 2 ) 


(Our specific choice of the binning will be given in Sec. 


IIC ) The error bars on r are then derived from the Fisher 


matrix via: 


K )- 2 



with "Ebb 1 = Cov [Cf which stands for the covari¬ 

ance matrix of the reconstructed, binned angular power 
spectrum of the B-mode. (Note that Cf denotes the 
estimator of the angular power spectrum, Cf.) 

The B-mode angular power spectrum as a function of 
r is modeled as: 


Cf (r)=rx7?+7^f, (4) 

with Tf and Tf&f two fiducial angular power spectra, 
which do not depend on r. The former is just obtained 
as the contribution of primordial gravity waves for r = 1 
(taking into account that the primordial B-mode is it¬ 
self lensed). The latter corresponds to the contribution 


of primary B-mode transferred into B-mode because of 
the gravitational lensing of large scale structures. We do 
not consider here a potential delensing of the B-mode 
anisotropies, and such a contribution will be assumed to 
act as an additional gaussian noise for the measurement 
of r. This is a simplifying assumption since the lensing- 
induced B-mode is non-gaussian leading to an additional, 
non-gaussian contribution to the covariance 551 . Gaus- 
sianity remains however a good approximation for band- 
powers which are narrow enough (A& < 100) [35], which 
is the case in our study. 

The covariance matrix Sf,b' is estimated using three 
different approaches, as described here. 


1. Mode-counting 


First, we rely on a naive mode-counting expression (or 
so-called / s k y -formula). In this case, the covariance on 
Cf is approximated by: 



2 5i,t ( c b , N e(fsk y ) 

(2£+l)f sky \ e + B\ 


(5) 

with N( the noise power spectrum, Be the beam of the 
telescope, and, f sky the portion of the celestial sphere, 
which is observed (or kept in the analysis). The noise 
power spectrum scales linearly with the sky coverage. 

The covariance matrix for the binned power spectrum 
is thus given by: 


Sbb' 


Y,(Pm) 2 * Cov (df,Cf) 

teb 


5b,b'- 


( 6 ) 


This is essentially used as a benchmark as such an eval¬ 
uation of the statistical error bars on the B-mode recon¬ 
struction underestimates the error bars coming from any 
numerical methods to be applied to the data. 


2. Minimum variance quadratic estimator 


Second we consider the error bars that could be in¬ 
curred by using a minimum variance quadratic estimator 
ESI HQl. The estimator is defined as follows: 


Cf 








(7) 

In the above, C = (dd^) is the covariance matrix of the 
maps of the Stokes parameter, and d is the column vector 
composed of (/, Q , U) (the trace operation is across 
pixels). The quantity N/j stands for the noise debias. 
Finally, F is the Fisher information matrix given by: 


F w 



r ac 3c 

dCf dCf 


( 8 ) 
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It is then shown that the covariance of the above esti¬ 
mator is given by the inverse of the Fisher matrix, i.e. 

Cov (cf,CPJ = F“^ ; . We remind that this estimator 

is precisely built to be the quadratic estimator with the 
lowest variance. 

If the 13-mode power spectrum is indeed estimated for 
each multipole, i (that is chosing A& = 1), this directly 
gives the following expression for the error bars expected 
on r: 


with the interpolation operator, Qbe- 

) 2tt 

W+ i) lf 1 e b ' (12) 

0 if l$b. 

The binned estimator, C^, is finally defined as: 

Cf = ^F-^,C> pt) (13) 


(<T r ) — F rr 



— C” 1 —c _1 

dr dr 


(9) 


with C the same covariance matrix but assuming that 
only Cf does depend on r, in line with our approach 
consisting in constraining the tensor-to-scalar ratio from 
the I3-mode’s measurements onljJ^J In the case of az- 
imutlrally symmetric patches, the numerical computation 
of such Fisher matrices (either F w or F rr ), can be per¬ 
formed in a reasonable time using the expression found 
in the appendix F of Ref. J32], and by using the S 2 hat 
package to perform spherical harmonic transforms 0U- 
IH]. (The use of this massively parallel package allows 
for a rapid computation of the covariance matrix for large 
sky coverages.). In the standard case, for brute force cal¬ 
culation the Fisher matrix requires 0{N^- ly ) operations 
to be evaluated, but in this calculation, evaluating the 
spin harmonics by recursion in i makes the computa¬ 
tional cost as 0(Ngm lnax ) 1 where Ng is the number of 
rings actually used. 

Practically speaking, one should nonetheless include 
the impact of binning, done as follows. First one defines 
the so-called optimal pseudo spectrum: 

Cj ° pt) = Tr 

One easily checks that ^(5^ opt ^ = 2 FwB^Cf, 

(where we also include the impact of an azimuthally sym¬ 
metric beam). One then introduces the matrix: 

PuFu’B^Qb'i'-, ( 11 ) 

ZGb Z'eb' 


d,,c - 1 S c - 


- Ne- (10) 


1 We note that from the complete definition of a r given by: 


(oy)" 2 = -Tr 


dCc-^c- 
dr dr 


with C the covariance matrix assuming that all the angular 
power spectra do depend on r, the equation is therefore re¬ 
placed by: 



In the above, the indices A, A' runs over TT, EE , BB , TE , TB 
and EB. The equation § is finally obtained assuming that only 
Cf in C does depend on the tensor-to-scalar ratio. 


with cT pt) = Yhz 1 Pvt' C'j° pt ' > the binned, optimal pseu¬ 
dospectrum. From that last definition, and making use 
of Eqs. ([7]) and ([8]), it is straightforward to show that: 

= [F-V [P Ml F v , P b , e ,] [(F- 1 )^,,,, (14) 

where summations over repeated indices (i.e. b \, b\ and 
l \, t\ ) is implicitly assumed, and f means the transpose 
operation. 

We note that this way of estimating the uncertainties 
on the power spectrum reconstruction is also relevant for 
maximum-likelihood approaches, see e.g. Ref. [ 35] . 


3. Pure pseudospectrum 

Third, we make use of the X 2 PURE code and Monte- 
Carlo simulations to estimate the covariance matrix ex¬ 
pected for the pure pseudospectrum approach. Details 
on the pure pseudospectrum estimator can be found in 
Refs. [221 [23] - In practice, the power spectrum is esti¬ 
mated within bandpower and the covariance matrix re¬ 
constructed from the MC simulations is directly S&f/. 
The numerical cost of this scales as O(N^). allowing 
for rapid MC simulations. 

For each sky coverage and for each value of r here- 
considered, we compute optimized sky apodizations to 
apply to the maps of Q and U. Those optimized sky 
apodizations are described in Refs. [231 El, and they 
allow for having the smallest error bars on the 13-mode 
power spectrum reconstruction within the context of pure 
pseudospectrum techniques. Those sky apodizations are 
a set of spin-0, spin-1 and spin-2 window functions to 
be applied to the maps of the Stokes parameter. They 
can be interpreted as the window functions, which make 
the pure pseudospectrum estimator as close as possible 
to the minimum-variance quadratic estimator [23]. 

Numerically speaking, computing those optimized sky 
apodizations may be long, especially for intricate shape 
of the observed region and/or low level of noise. Using 
an iterative method, the numerical cost is C>(WterlV 2 ix ), 
with IVjter a number of iterations ranging from few tens to 
few hundreds for simple patch geometry and noise level 
considered here (see Sec. Ill in [24] )■ We stress that for 
a given sky patch, those sky apodizations are to be opti¬ 
mized for each value of the tensor-to-scalar ratio and bin- 
by-bin. Taking into account the number of bins (see Sec. 
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, the number of sky fractions and the number of val¬ 
ues of r, which are sampled in this study, this means that 
1664 of such sky apodizations have to be computed. For¬ 
tunately, in the case of homogeneous noise and patches 
with relatively simple contours (which is obviously the 
case for a spherical cap), it was demonstrated in Ref. 
}24) that an approximated but numerically fast compu¬ 
tation of those sky apodizations in the harmonic domain 
is possible, and indeed leads to error bars equal to those 
obtained thanks to a direct, pixel-based computation of 
the optimized sky apodizations. The numerical cost of 
this technique is reduced to O(N^) which allows us to 
derive optimized sky apodizations for each value of the 
sky coverage, and for each value of the tensor-to-scalar 
ratio. 


C. Power spectrum uncertainties 

Our chosen bandpowers for reconstructing Cf are the 
following. 

The first bin starts at £ = 20 and we use a con¬ 
stant bandwidth, A ft = 40. Our last bin extends up to 
£ = 1020. The value of the maximum multipole is chosen 
in order to include all the relevant contributions from the 
primordial 13-mode, that is until the lensing 13-mode be 
the dominant contribution to the total 13-mode power. In 
addition, for the experimental cases under consideration 
in this paper, we use beamwidths up to 8 arcminutes, 
corresponding to a cut-off of £ ~ 1300. The choice of the 
bandwidth is mainly motivated by the use of the pure 
pseudospectrum estimator. Especially, it is mandatory 
for the numerical inversion of the mode-mixing matrices 
to be possible. We also note that the bandwidth is wide 
enough so that the correlations between different bins 
are nearly uncorrelated in the covariance of pure pseu¬ 
dospectrum estimator. We discuss the role of the bin 
width later on. We stress that the multipoles ranging 
from £ = 2 to £ = 20 (corresponding to the reionization 
peak, and gathered in one bandpower) are actually used 
in the pure pseudo-CV estimation of Cf. However, given 
the limited sky coverages considered here such low mul¬ 
tipoles are difficult to estimate and hardly constrained 
by the data. This bin is therefore not included in our 
analysis of the signal-to-noise ratio. 

The uncertainties on the estimated power spectrum of 
the 13-mode as functions of the sky fraction, are shown 
in Fig. [T] where four selected values of the sky cover¬ 
age are depicted: 1%, 3.5%, 7% and 10%. The tensor- 
to-scalar ratio chosen for this figure is r = 0.1. Each 
panel corresponds to a different approach to derive the 
covariance matrix, 'Sibb'- mode-counting, minimum vari¬ 
ance quadratic estimator, and, pure pseudospectrum es¬ 
timator (from left to right). 

As expected, we do observe that the lowest error bars 
are the ones from the mode-counting estimation of the 
uncertainties, while the highest error bars are obtained 
from the pure pseudo-CV estimator. The error bars from 


the minimum variance, quadratic estimator lie between 
those two. At the largest accessible scales, 2 < £ < 100, 
the error bars from the pure pseudospectrum estimator 
are ~ 1.25 greater than the optimistic mode-counting 
estimation. Similarly, the error bars from the pure 
pseudospectrum estimators are at most ~ 1.1 higher 
than the ones derived from the minimum variance, 
quadratic estimators. At the smaller angular scales 
where lensing dominates, the three approaches lead to 
almost the same uncertainties. 

The behaviour of the uncertainties as a function of the 
sky fraction is common to the three approaches. At the 
smaller angular scales first (for multipoles greater than 
~ 100), the behavior is monotonic since the uncertainties 
systematically increase with the value of f s k y . This is be¬ 
cause at these scales, the variance is dominated by the 
noise, which increases with the sky fraction. At larger 
scales however (for multipoles smaller than ~ 100), the 
uncertainties have a more intricate behaviour. First one 
notes that the uncertainties decrease from / s k y = 1% to 
/sky = 3.5%. This is because the variance is dominated 
by sampling variance, which decreases for higher values 
of /sky- Second, one notes that uncertainties at £ < 100 
then increases for a sky coverage ranging from 3.5% to 
10%. This means that for f s k y > 3.5%, the noise is 
now dominating the variance. Once this transition value 
of / s ky ~ 3.5% is crossed, the noise contribution domi¬ 
nates the variance for all our considered angular scales, 
£ € [20,1020]. Therefore, the variance will monotoni- 
cally increase with / s k y at all the relevant angular scales 
once / s ky > 3.5%. (Note that an identical behaviour is 
observed for the other values of r, though the specific 
value of /sky at which the transition occurs depends on 
the specific value of r.) 


D. Signal-to-noise ratio on r 

1. Numerical results 

The signal-to-noise ratio on r is computing using Eq. 
© considering the three above-described methods to es¬ 
timate the uncertainties on the i3-mode reconstruction, 
"Ebb 1 ■ We remind that the summation in ^ is performed 
over bandpowers with a bandwidth of Af, = 40 and con¬ 
sidering a range of multipoles from £ = 20 to £ = 1020. 

Our numerical results on the signal-to-noise ratio for r 
are gathered in Fig. [2j showing (S/N) r as a function of 
the sky coverage. Each panel corresponds to a given value 
of the tensor-to-scalar ratio, r = 0.07, 0.1, 0.15, and, 0.2 
(from top to bottom). For each panel, the black, red, 
and blue crosses correspond to the signal-to-noise ratio 
derived by using the mode-counting, the minimum vari¬ 
ance quadratic estimator, and, the pure pseudo-C^ esti¬ 
mator, respectively. The horizontal, dashed line marks a 
3cr detection. The sky fraction varies from 0.1% to 10%, 
what is wide enough to sample the maximal values of the 
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FIG. 1. Uncertainties on the reconstructed B-mode power spectrum (dashed color curves) for four values of /sky: 1%, 3.5%, 
7% and 10%. The solid black curve is the input Cf for r = 0.1. Each panel corresponds to a different approach to derive the 
covariance matrix £ w. mode-counting, minimum variance quadratic estimator, and, pure pseudospectrum estimator (from 
left to right). 


signal-to-noise ratio. We note that the signal-to-noise 
ratio keeps decreasing for / sky > 10%. This is because 
for the level of noise and values of r here considered, 
the uncertainties on the reconstructed B-mode are noise 
dominated at all scales for / sky > 10%. Similarly, the 
(S/N) r keeps decreasing for / sky < 0.5%, because the 
uncertainties on angular scales greater than a degree are 
dominated by the sampling variance for such low values 
of the sky fraction. 

For the case of the mode-counting first, the signal-to- 
noise ratio is systematically greater than 3 for all the 
considered values of the sky coverage, and for all the 
considered values of the tensor-to-scalar ratio. 

Considering then the case of the minimum variance, 
quadratic estimator, the signal-to-noise ratio on r is sys¬ 
tematically greater than 3 for r = 0.1, 0.15 and 0.2. For 
a tensor-to-scalr ratio of r = 0.07, the (S/N) r is greater 
or equal to three for 1% < / sky < 9%. 

Assuming finally a pure pseudospectrum recon¬ 
struction of the B-mode, the signal-to-noise ratio is 
systematically greater than 3 for r = 0.15 and r = 0.2 
only. For a tensor-to-scalar ratio of r = 0.1, a mea¬ 
surement of it with a (S/N) r of at least 3, is possible 
for / sky > 1%. For a smaller value of r = 0.07, its 
measurement with (S/N) r > 3 is possible assuming 
1.5% < / sky < 7%. We note however here that for 
r = 0.07 and r = 0.1, the signal-to-noise ratios remains 
greater than 2 assuming a pure pseudospectrum recon¬ 
struction of Cf. 

As expected from the error bars on the reconstructed 
B-mode, the highest and lowest (S/N) r ’s are respectively 
obtained from the mode-counting estimation, and the 
pure pseudospectrum estimator, while the (S/N) r from 
the minimum variance quadratic estimator lies between 
those two. This is the case for all the values of the tensor- 
to-scalar ratio we consider. At the peak, the signal-to- 


noise ratio from the pure pseudospectrum estimation of 
the Cf is ~15% (r = 0.2) to ^20% (r = 0.07) smaller 
than the signal-to-noise ratio derived from the optimistic 
mode-counting. This means that the statistical signifi¬ 
cance on the measurement of r by using the optimistic 
mode-counting is overestimated by a factor ~ 1.25 as 
compared to the more realistic case of the pure pseu¬ 
doreconstruction of the B-rnode. 

Similarly, the (S/N) r from the pure pseudospectrum 
estimator is ~1.5% (r = 0.2) to ~8% (r = 0.07) smaller 
than the signal-to-noise ratio derived from the minimum 
variance, quadratic estimator. Using the minimum vari¬ 
ance, quadratic estimator to estimate the B-rnode, as 
compared to the use of the pure pseudospectrum, thus 
translates into a gain in the statistical significance on 
the measurement of r, of a factor 1.01 to 1.08. This gain 
appears rather small but is larger for smaller values of 
the tensor-to-scalar ratio. 


2. Optimization of the sky coverage 

As clearly shown in Figs. [2j there exists a value of the 
sky coverage, which maximizes the signal-to-noise ratio 
on r. This optimal value of / sky was already observed 
in Ref. [IS], using only the mode-counting expression 
for the statistical error bars on the B-rnode estimation 
though. We found that such an optimal value also exists 
using the minimum variance quadratic estimator or the 
pure pseudo-C^ estimator. This is intuitively understood 
as follows. The statistical uncertainties on the angular 
power spectrum estimation have two sources, the sam¬ 
pling variance, which is dominant at the largest angular 
scales, and the noise variance dominating at the smallest 
angular scales. Reducing the sampling variance is ob¬ 
tained by covering a large fraction of the sky. However, 
for a given sensitivity and a given time of observation, 
covering a large fraction of the sky inevitably translates 





7 







+ ++++ +H 

h +++. 



X 1 

+ + 
+ 

+ 

h + + + H 
+ ^ 
H + + M 

E*tt< 

k '""+""h 

Lj._h 

h + + 

+ + 

+ 



+ r4 
+ + Mod 
+ + QML 

-. 07 + X 
e counting 




+ + Pur 



0 2 4 6 8 10 


Fraction of sky in % 









H 

+ 

h + + + H 

H'+' + H 

h + 1 

H + H 

h 

+ H 

+ t< 



*: i 

H 

h . 

+ + 

+ 






+ 



r=0.15 




+ + out 
+ + Pur 

e method 




6 8 10 


Fraction of sky in % 


9- 

8 

H 

7 + 


H + + + H 

' + 1 

h + H 

H 

+ H 

6 + j 


^ + 1 

t t i 

1 = 

5+ ++ 






+ + 






4 + 



r=C 
+ + Moc 

.2 

e countin 

g 




+ + QM 

+ + Pur 

e method 






io 


Fraction of sky in % 


FIG. 2. Signal-to-noise ratio for the estimation of r from B- 
mode polarization data shown as a function of the sky cover¬ 
age. The uncertainties are computed using three different ap¬ 
proaches: mode-counting (black crosses), minimum variance 
quadratic estimator (red crosses) and pure pseudospectrum 
estimator (blue crosses). Each panel corresponds to a different 
fiducial value of the tensor-to-scalar ratio, r = 0.07, 0.1, 0.15 
and 0.2 from top to bottom. 


r 

0.07 

0.1 

0.15 

0.2 

/IT 3 [%] = 

Mode-counting 

2.0 

3.0 

4.0 

5.0 

Minimum-variance Ct 

2.5 

2.5 

3.5 

5.0 

Pure pseudo-CV 

3.5 

3.5 

5.0 

4.0 

(S/N)(. opt) : 





Mode-counting 

4.4 

5.5 

7.0 

8.2 

Minimum-variance Ci 

3.7 

4.7 

5.9 

7.0 

Pure pseudo-G^ 

3.4 

4.4 

5.8 

6.9 


TABLE I. Values of the sky fraction maximizing the signal- 
to-noise ratio on r, and maximal values of the signal- 

to-noise ratio, (S/N)r° pt ' ) . This is given for each techniques 
used to estimate the uncertainties on the reconstruction of 
the B-mode. We notice that concerning the pure method, 
the position of the maximum value is not as defined as the 
other methods (see text). 


into a higher level of noise per pixel. One should therefore 
find the good balance between sampling and noise vari¬ 
ance so as to minimize the total error on given targetted 
parameters, which is r here. 

The salient features of those results are summarized in 
the table HI For each value of the tensor-to-scalar ratio 
and for each techniques used to compute uncertainties 
on the 73-mode, we provide the values of the sky fraction 
maximizing the signal-to-noise ratio, / s ^ pt \ Its associ¬ 
ated (thus maximal) value of the signal-to-noise ratio, 
(S/N)[ opt ^ is also reported in this table. We stress that 
the position of the peak of (S/N) r is well defined for the 
mode-counting and the minimum variance, quadratic es¬ 
timator. Such a position of the peaking signal-to-noise 
ratio is however less pronounced for the case of the pure 
pseudo-CV estimation of the 73-mode (see e.g. the case 
r = 0.7 for which a range of 2% < / s k y < 6% leads 
roughly to the same (S/N) r ). This means that the values 
of f^y T> reported in Tab. [i] for the case of the pure pseu¬ 
dospectrum approach are more indicative than a sharply 
defined value. 

For all the approaches used to estimate the uncertain¬ 
ties on the 73-mode, we observe that the optimal sky 
fraction increases with the value of the tensor-to-scalar 
ratio. This is because for higher values of r, the signal 
in the 73-mode is higher. One should therefore minimize 
first the sampling variance by increasing the observed 
part of the sky. 

Except for the case of r = 0.2, we note that the optimal 
sky coverage assuming a minimum variance, quadratic es¬ 
timator slightly differs by 0.5% (either higher or lower) 
than the value of / s ^ pt ' ) (in %) as inferred from the mode¬ 
counting. We also note that the optimal sky fraction 
obtained for the pure pseudo-reconstruction of the 73- 
mode differs by 1% to 1.5% (depending on the value of r) 
from the one inferred from the mode-counting estimation 




































of the uncertainties on the £?-mode. Nevertheless, the 
values of the sky fraction for which the detection of the 
tensor-to-scalar ratio is peaking in the case of the mode¬ 
counting expression and the minimum-variance quadratic 
estimator fall in the range of optimized / s k y as derived 
from the pure pseudo-C^ estimator. Those numerical re¬ 
sults therefore show that (at least) for the range of values 
of r here-considered, the value of the sky coverage, which 
maximizes the measurement of the tensor-to-scalar ratio 
is rather independent on the adopted method for evalu¬ 
ating the statistical uncertainties on the 33-mode recon¬ 
struction. This means that using the mode-counting ex¬ 
pression, though underestimating the error bars, allows 
for a rapid and reliable search of the range of values of the 
optimized sky fraction. (Obviously, such an optimization 
of / s ky based on the mode-counting expression is reliable 
providing the final data set to be analyzed using either 
the minimum-variance quadratic estimator or the pure 
pseudospectrum estimator.) 


3. Impact of binning 


For the two specific cases of the mode-counting uncer¬ 
tainties and the minimum variance, quadratic estimator, 
we note that an explicit reconstruction of the power spec¬ 
trum is not mandatory to derive the (S/N) r in the Fisher 
formalism. One can indeed directly plugged in Eq. ([3]) 
the formulas (f6l) or (14). This allows for a study of the 


impact of binning on the signal-to-noise ratio, letting the 
bandwidth to vary from A 6 = 1 (i.e. no binning) to 
A b = 40 (i.e. the binning imposed by the use of the 
pseudospectrum estimator in this analysis). 

The impact of binning is illustrated in Fig. [3] show¬ 
ing the signal-to-noise ratio on r = 0.1 as a function 
of the sky coverage. The grey (red) area corresponds 
to the (S/N) r using the mode-counting (minimum vari¬ 
ance quadratic estimator) to estimate the uncertainties 
the angular power spectrum of the 33-mode. For each 
shaded area, the highest signal-to-noise ratio is obtained 
for A 5 = 1 and the lowest for Aj, = 40. As a reference, we 
also show the (S/N) r obtained with pure pseudospectrum 
reconstruction (thus using a bandwidth of Aj = 40) de¬ 
picted by the blue crosses. The overall effect of increasing 
the width of the bandpower is to lower the signal-to-noise 
ratio. The decrease is however more pronounced for the 
case of the minimum variance, quadratic estimator than 
for the mode-counting estimation of the error bars on 
the reconstructed 13-mode. This is due to the fact that 
correlations between multipoles (or bandpowers) are ac¬ 
counted for in the minimum variance, quadratic estima¬ 
tor, while those are supposed to be systematically van¬ 
ishing for the mode-counting estimation of the covariance 
matrix. This additional piece of information contained 
in the correlations is therefore partially lost by averaging 
over bandpowers. We also checked that artificially impos¬ 
ing those off-diagonal correlations to be zero lowered the 
signal-to-noise ratio in the minimum variance method, 
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Fraction of sky in 


-r=Q.T-. 

- - Mode counting 
-- QML 

+ v Pure method 


FIG. 3. Signal-to-noise ratio for r = 0.1 as a function of 
the observed sky fraction, derived for three methods used 
to estimate the uncertainties on the 33-mode reconstruction: 
mode-counting (black area), minimum variance quadratic es¬ 
timator (red area), and, pure pseudospectrum estimator (blue 
crosses). For the two first methods, we let the bandwidth of 
the bins to vary from Aj, = 1 (highest (S/N) r ) to Aj, = 40 
(lowest (S/N) r ). For the specific case of the pure pseudospec¬ 
trum estimator, the reconstruction of the Cf requires to 
use the bandwidth A& = 40. (We remind that the range 
of multipoles used to compute the signal-to-noise ratio is 
20 < I < 1020 .) 


although we note that once the bins are sufficiently wide 
the effect of the binwidth on the (S/N) r should be weak. 

The maximum values of the (S/N) r obtained using a 
bandwidth of A;, = 1, and a bandwidth of Af, = 40, 
are reported in Tab. [nj for each values of r and for the 
mode-counting and the minimum variance quadratic es¬ 
timator. For each cases, we also report the value of the 
sky coverage corresponding to that maximum. For the 
mode-counting approach, increasing the bandwidth from 
Aj = 1 to A 6 = 40, degrades the maximum (S/N) r by 
a factor ~ 9% for r = 0.07 and 0.1, and, by a factor 
~ 2% for r = 0.15 and 0.2. This however only mildly 
affects the values of the sky fraction at which the maxi¬ 
mum is achieved. The impact of binning is more marked 
for the minimum variance quadratic estimator however. 
Increasing the bandwidth from Aj, = 1 to Aj = 40, here 
degrades the maximum (S/N) r by a factor ~ 11% for all 
the values of the tensor-to-scalar ratio considered in this 
study. Similarly, the values of the sky fraction (in %) at 
which this maximum is achieved is systematically low¬ 
ered (except for the case r = 0.7), by 1% for r = 0.1 and 
by 2% for r = 0.2. We note that despite these changes 
in the value of with the bandwidth, the optimized 

values of the sky fraction still fall in the range of opti¬ 
mized / s ky as derived from the pure pseudo- 6 / estimator. 
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r 

0.07 

0.1 

0.15 

0.2 

(S/N)j° pt) and / s ( ° y pt) : 
mode-counting: 

= 1 

Af, = 40 

4.6 (2.5%) 
4.4 (2.0%) 

5.7 (3%) 

5.5 (3%) 

7.2 (4%) 

7.0 (4%) 

8.4 (5.5%) 

8.2 (5%) 

minimum-variance Cf. 

At = 1 

A,, = 40 

4.2 (2.5%) 
3.7 (2.5%) 

5.3 (3.5%) 
4.7 (2.5%) 

6.8 (5%) 

5.9 (3.5%) 

8.0 (7%) 

7.0 (5%) 


TABLE II. Maximum values of the (S/N) r for the mode-counting method (upper part), and the minimum variance quadratic 
estimator (lower part). This is derived assuming no-binning (i.e. At = 1), or using a binning with a bandwidth of At = 40. In 
parentheses is the values of / s ky (in percent) at which this maximum is reached. 


III. MEASURING THE TENSOR-TO-SCALAR 
RATIO: SELECTED EXAMPLES 

A. Experimental setups 

We turn to the question of the detection of r in more 
realistic cases. Clearly, a spherical cap is ideal. The is¬ 
sue of leakages is strongly related to the detailed shape 
of the contours of the observed (or kept-in-the-analysis) 
portion of the sky (see e.g. the figure 20 of Ref. [24] for 
the impact of the shape of the mask on the statistical 
error bars). A spherical cap then leads to the smallest 
amount of leakages for a given sky fraction since its con¬ 
tour has the smallest perimeter for that given sky frac¬ 
tion. To this end, we consider three archetypal cases, 
which capture the main characteristics of ongoing, or 
being-deployed, small-scale experiments (ground-based 
or balloon-borne), a possible upgrade of those ground- 
based experiments to an array covering a rather large 
fraction of the sky (~ 50%), and, a possible satellite mis¬ 
sion covering the entire celestial sphere. 

First, we consider a part of the sky, which is typi¬ 
cal of small-scale experiments, inspired by the design 
study of the ebex experiment [47] and very similar to the 
patch and noise level used in Ref. [4Sj mimicking current 
ground-based efforts such as ACTpol or polarbear. 
This roughly covers 1% of the sky with a highly inho¬ 
mogeneous noise, as shown in the left panel of Fig. |4] 
The average noise level is ~ 5.75/rK-arcminute and the 
beamwidth is 8 arcminutes. This setup will be referred 
as ’’small-scale experiment” hereafter. 

Second, we consider the potential case of an array 
of ground-based telescopes allowing for covering ~ 50% 
of the celestial sphere with a (here-considered homoge¬ 
neous) noise of 1/rK-arcminute and a beamwidth of 3 
arcminutes. These specifications roughly correspond to 
the Stage IV of future CMB experiment as reported in 
e.g. Ref. [32]. For such a large fraction of the sky, mask¬ 
ing the regions with high foreground galactic emissions 
is needed. To this end, we consider that the entire galac¬ 
tic south hemisphere would be observed and we apply a 
galactic mask, and a mask for point sources. We make 


use of the R9 galactic mask used for polarized data of 
WMAP and add the point-sources mask [50]. The result¬ 
ing sky coverage is depicted in the middle panel of Fig. [4j 
and it roughly covers ~36% of the sky. (We note that this 
sky coverage is just the restriction to the galactic south 
hemisphere of the kept-in-the-analysis portion of the sky 
for a satellite mission as defined below.) This config¬ 
uration will be referred as ’’array of telescopes” in the 
following. One could have chosen also the south ecliptic 
hemisphere, however while this consideration is impor¬ 
tant for e.g. scanning strategy and foreground treatment, 
this is less relevant for the minimization of the E-to-B 
leakage for which one key point is the complexity of the 
contours of the patch, which would be reflected in both 
choices of observation. 

Third, we consider a potential satellite-like mission 
with homogeneous noise at a level of 2.2p,K-arcminute 
and a beamwidth of 8 arcminutes m- The portion of 
the sky to be analysed is displayed in the right panel of 
Fig. @ It corresponds to 71% of the entire celestial 
sphere. As for the previous case, the removed portion 
of the sky corresponds to the WMAP galactic and point- 
sources mask. This will be referred as ’’satellite mission” 
hereafter. 

B. Sky apodizations 

For each of those cases, the uncertainties on the re¬ 
constructed angular power spectrum are obtained using 
the pure pseudospectrum estimator, using the same set 
of bandpowers as in Sec. |II C| We note however that the 
first bandpower, ranging from £ = 2 to £ = 20 will be 
explictly shown hereafter since those scales are now ac¬ 
cessible for the case of an array of telescopes, and the 
case of a satellite mission. 

Optimizing the sky apodization to be applied to the 
maps is a key step to reach better performances on the 
reconstructed Cf. Two classes of sky apodizations have 
been proposed in the literature (see Refs. [231 !23])- The 
first class consists in an analytical formula fullfilling 
the appropriate Neuman and Dirichlet boundary condi¬ 
tions using arches of a sine function. The apodization 
length can be subsequently optimized, bin-per-bin, to 
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Cut distribution of sky observations 



FIG. 4. Left Panel: A potential noise distribution (number of hits per pixel) for small-scale experiments covering ~1% of the 
sky. The noisiest pixels have been discarded and the noise distribution ranges from 10 5 to 10'. Middle panel: Portion of the 
sky kept for estimating the B-mode angular power spectra for a possible array of ground-based telescopes. It roughly covers 
36% of the sky (Note that for this analysis the choice of galactic or ecliptic south hemisphere is equivalent, see text for more 
information). Right panel: Portion of the sky kept for estimating the B-mode angular power spectra for a possible satellite 
mission, roughly covering 71% of the sky. 


minimize the uncertainties on the estimated Cf. The 
second class is a set of minimum-variance optimized sky 
apodizations. They can be computed either in the pixel 
domain as originally proposed in S5] (thus allowing for 
relaxing the boundary conditions and keeping track of 
the informations about the B-modes contained in the 
so-called ambiguous modes), or in the harmonic domain 
as proposed in [24] providing the noise is homogeneous 
but here forcing the Neuman and Dirichlet boundary 
conditions to hold. The pixel-based computation is more 
flexible and general (being e.g. applicable to cases of 
inhomogeneous noise), and it was shown to generically 
lead to better performances. It is however numerically 
costly, while the harmonic-based computation is very 
rapiqj 

To illustrate the specific case of the minimum-variance 
sky apodizations (computed in the pixel domain), an ex¬ 
ample of its scalar, spin-0 component is depicted in Fig. 
[5]for the three experimental setups. This sky apodization 
has been optimized for the bandpower comprised between 
£ = 60 and £ = 100 and for a value of the tensor-to- 
scalar ratio r = 0.05. We note that the sky apodization 
computed on the sole galactic south hemisphere for an 
array of telescopes (see middle panel of Fig. [5]) is not 
the restriction to the galactic south hemisphere of the 
sky apodization for a satellite mission computed on both 
hemispheres (see right panel of Fig. [5|. (This would ob¬ 
viously be the case for analytic sky apodizations.) This 
is mainly due to the fact that the global shape of the 
observed portion of the sky is taken into account in the 


2 This is because those optimized sky apodization are obtained as 
the solution of linear system of size JV p ; x . This linear system cor¬ 
responds to a convolution, which is drastically simplified in the 
harmonic domain (the convolution kernels becoming diagonal) 
assuming the noise is homogeneous. 


optimization process. This is clearly seen by noticing 
that for a typical patch covering both hemispheres, there 
is a rather large apodization length parallel to the az¬ 
imuthal direction while the apodization length along the 
zenithal direction is rather small. On the contrary for 
a typical patch covering only the south hemisphere, the 
apodization length is essentially along the zenithal di¬ 
rection. Such a difference is explained by the important 
hole along the zenithal direction in the north hemisphere, 
which leaves its footprint on the final shape of the sky 
apodization. (There is a second source of differences due 
to the different beamwidth and noise level, which are 
considered for both types of experiments. They however 
have a tiny impact on the orientation of the apodization 
length since they are evenly distributed over the sky and 
therefore do not have any preferred directions.) 

A couple of comments about the numerical computa¬ 
tion of the pixel-based, minimum-variance sky apodiza¬ 
tions are in order here. They are theoretically built to 
give the smallest uncertainties in the context of the pure 
pseudospectrum estimators. However, they are pratically 
computed from a Preconditionned Conjugate Gradient 
(PCG) algorithm, which efficiency strongly depends on 
the experimental configurations, especially with respect 
to the noise level and its distribution over the patch, as 
well as with respect to the complexity of the contours of 
the mask. 

First, the number of iterations in our implemented 
PCG rapidly increases for lower levels of noise: at the 
largest angular scales (£ < 20) the number of iterations 
ranges from ~ 100 for a noise level of 5.75p,K-arcminute 
to ~ 3000 for a noise level of l^K-arminute (the num¬ 
ber of iterations being one order of magnitude smaller 
for smaller angular scales, £ > 20). The B-mode angular 
power spectrum is estimated using the same binning as 
in the previous section, leading to N^ ln = 26, and we se¬ 
lected 6 values of r. Considering 3 experimental setups, 
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Spin = 0 Window 


FIG. 5. Spin-0 (scalar) component of the optimized sky apodization for a small-scale experiment covering ~1% of the sky with 
inhomogeneous noise (left panel), for a possible array of ground-based telescope covering ~36% of the sky with homogeneous 
noise (middle panel), and, for a possible satellite mission covering ~71% of the sky with homogeneous noise (right panel). This 
sky apodization is optimized for a bandpower ranging from i = 60 to l = 100 and for a value of the tensor-to-scalar ratio 
r = 0.05. 


this would translate into ~500 optimized sky apodization 
to compute, which is numerically too costly. We there¬ 
fore compute the optimized sky apodization for r = 0.05 
only but use them for all the here-considered values of 
r, meaning that the signal-to-noise ratios obtained for 
r ^ 0.05 may be suboptimal within the context of the 
pure pseudospectrum estimatoi]^] 

Second, it is not guaranteed that the algorithm con¬ 
verges towards the optimal solution, especially for inho¬ 
mogeneous noise (see Sec. IV C of Ref. [24] where it 
was shown that trimming out the external, noisiest pix¬ 
els is required), or a very low level of noise (see Ref. [23] 
mentionning that convergence is not reached for a noise 
level of ~ 1/rK-arminute, corresponding to the level of 
the array of telescopes case). This means that the per¬ 
formances of those sky apodizations have to be assessed 
using numerical simulations at the level of power spec¬ 
trum reconstruction, comparing the resulting error bars 
on the estimated power spectra to the error bars obtained 
by using the other types of sky apodizations. 

C. Power spectrum uncertainties 

The relative performances of the different sky apodiza¬ 
tions are appraised at the level of power spectrum un¬ 
certainties. For each case we performed a series of 500 
Monte-Carlo simulations to compute the statistical un¬ 
certainties on the reconstructed R-mode angular power 
spectra, assuming different kinds of sky apodizations. 


3 In the specific case of small-scale experiments, it was however 
shown in Ref. m that the optimization process is mainly driven 
by the noise level and the amount of .E-modes leaking into the 
E-mode, and poorly affected by the amplitude of the primordial 
E-mode. This means that the resulting sky apodizations may be 
mildly dependant on the assumed value of r and that the derived 
signal-to-noise ratios are only slightly suboptimal for the case of 
small-scale experiments. 


Such performances have been exhaustively studied for the 
small-scale experiment case and the satellite mission case 
(see Refs. [23] [25 and Ref. [30], resp.). On the contrary, 
the applicability of the pure pseudospectrum estimator 
for the case of an array of telescopes was hitherto not 
studied. We then performed numerical simulations using 
the different classes of sky apodizations to assess the effi¬ 
ciency of the pure pseudospectrum reconstruction of the 
R-mode, and subsequently select those sky apodizations, 
which lead to the smallest uncertainties. 

In this section, we only briefly review the major con¬ 
clusions concerning the cases of a small-scale experiment 
and a satellite mission. Then, we present the results of 
our numerical investigations for the case of an array of 
telescopes. 

1. Optimal apodizations: small-scale experiments and 
satellites missions 

For the small-scale experiment case, it was shown 
that the lowest uncertainties in the range £ £ [2,1020] 
were obtained using either the pixel-based optimized sky 
apodizations or analytic sky apodizations appropriately 
chosen to minimize the variance per each bin. However, 
the harmonic-based computation of the sky apodization 
fails in providing error bars comparable to the previous 
ones in the entire range of multipoles considered here, 
simply because the noise is inhomogeneous. (We refer to 
Fig. 24 of Ref. [25 and discussions therein.) 

For the satellite mission case, the pixel-based computa¬ 
tion of the minimum-variance sky apodization yields the 
smallest uncertainties for the range i £ [2,1020]. Simi¬ 
lar performance is obtained by using the harmonic-based 
computation of these sky apodizations for 150 < i < 600. 
The error bars however drastically increase for larger an¬ 
gular scales, the reason for that being the intricate con¬ 
tours of the galactic mask and the point-sources mask, 
which require to relax the Neuman and Diriclilet bound¬ 
ary conditions to keep (part of) the informations about 
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the B-mode contained in the ambiguous modes. For such 
a case, analytic sky apodizations fails in providing com¬ 
parable error bars at the largest angular scales. (We refer 
to Fig. 4 of Ref. [30] and discussions therein.) 

2. Optimal apodizations: arrays of telescopes 

For the case of an array of telescopes, we systemati¬ 
cally search for the type of sky apodizations, which lead 
to the smallest uncertainties bin per bin and for each val¬ 
ues of r considered in this study. We first found that the 
harmonic-based, optimized sky apodizations yield to er¬ 
ror bars higher than the analytic sky apodizations or the 
pixel-based, optimized sky apodizations. This is similar 
to what was observed for the case of a satellite mission 
and the inefficiency of the harmonic-based, optimized sky 
apodizations is due to the complexity of the contours of 
the mask. 

An example of the uncertainties for r = 0.1 and 
using the pixel-based, optimized sky apodizations (the 
dashed-red curve) or the analytic ones for different values 
of the apodization length (the dashed-blue curves) is 
shown in Fig. [6] The solid-black curve stands for the 
input angular power spectrum and the dashed-black 
curve stands for the binned, mode-counting computation 
of the error bars. This first shows that the pixel-based, 
optimized sky apodizations perform the best at the 
largest angular scales. This is systematically so for the 
two first bins. For the third bin, the pixel-based, sky 
apodizations and the analytic ones perform the same for 
r > 0.1, while for r < 0.1, the analytic sky apodizations 
with an apodization length of 4 degrees lead to a 
smaller error bar than the pixel-based, optimized sky 
apodization. However at smaller scales, £ > 100, the 
smallest error bars are systematically obtained by using 
an analytic sky apodizations with an apodization length 
of 2 degrees for intermediate scales, 100 < £ < 300, 
and an apodization length of 1 degree for small 
scales, i > 300. We found this to be independant of the 
value of r (at least for the grid of values considered here). 

We note that the apparent failure of the pixel-based, 
optimized sky apodizations here may be rather due to 
practical difficulties in computing such apodizations suf¬ 
ficiently accurately, rather than an indication of some 
fundamental problems. Indeed, we have found that for 
the noise levels the iterative solver used to compute 
the apodizations converges extremely slowly (as also ob¬ 
served in Ref. 23]) potentially preventing us in practice 
from achieving sufficient precision. 


3. Summary on the power spectrum uncertainties 

As a summary, the smallest statistical uncertainties 
obtained for r = 0.1 are shown in Fig. [7] in which the or¬ 
ange, red and burgundy curves stand for the small-scale 
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FIG. 6. Statistical uncertainties on the reconstructed B-mode 
for the case of an array of telescopes. The solid-black curve 
stands for the input angular power spectrum with a tensor- 
to-scalar ratio equal to r = 0.1. The dashed-black curve is 
for the binned, mode-counting uncertainties used as a bench¬ 
mark. The error bars obtained by using the pixel-based, op¬ 
timized sky apodizations (called PCG in the figure legend) is 
represented by the dashed-red curve. The different dashed- 
blue curves correspond to the error bars obtained by using 
analytic sky apodizations with an apodization length ranging 
from 1 degree to 6 degrees. 


experiment, an array of telescopes and a satellite mission 
respectively. For each experimental setup, we show the 
smallest error bars, which are attained for each band- 
power. For the cases of a small-scale experiment and 
a satellite mission, this is obtained by using the pixel- 
based, optimized sky apodizations throughout the entire 
range of angular scales. For the case of an array of tele¬ 
scopes, the pixel-based, optimized sky apodizations are 
used for multipoles lower than 100, while analytic sky 
apodizations with an apodization of 2 degrees and 1 de¬ 
gree are used in the range 100 < £ < 300 and in the range 
i > 300, respectively. 

As expected, the higher uncertainties are the ones from 
a small-scale experiment due to the tiny fraction of the 
sky it covers, and the relatively high level of instrumental 
noise. We provide the uncertainties for the first band- 
power, 2 < £ < 20, for completeness. These scales are 
nonetheless unaccessible starting from a map covering 1% 
of the sky due to the high uncertainties, as already stated 
in Sec. [Hi 

For angular scales going from £ = 2 to £ ~ 600, the 
smallest error bars corresponds to a satellite mission. 
This is because at these angular scales, the uncertain¬ 
ties are dominated by sampling variance and a satellite 
mission, as compared to an array of telescopes, benefits 
from its larger sky fraction. For multipoles smaller than 
~ 100, the error bars from a satellite mission are roughly 
1.5 to 2 times smaller than the error bars obtained from 
an array of telescopes, which is in line with the fact that 
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FIG. 7. Statistical uncertainties on the reconstructed 13- 
mode’s angular power spectra with the pure pseudospectrum 
estimator. The orange, red and burgundy lines stand for the 
small-scale experiment, an array of telescopes and a satellite 
mission respectively. The black line corresponds to the input 
angular power spectrum for r = 0.1. 


the sky fraction for a satellite mission is 2 times higher 
than the sky fraction observed by an array of telescopes, 
thus reducing the error bars by a factor of ~ y/2 as com¬ 
pared to the error bars from an array of telescope. 

Nevertheless, at small scales, l > 600, smaller error 
bars are obtained from an observation by an array of tele¬ 
scopes. This is because in that regime, the uncertainties 
for the case of a satellite mission are dominated by the 
noise term, Nt/Bj. Since the noise for a satellite mission 
is four times higher (in power spectrum) than the noise 
for an array of telescopes, and the beam is more than 
two times higher, this increase of the variance largely 
overcome the decrease due to a larger sky coverage. This 
quantitatively explains why at those small angular scales, 
the lowest error bars on the 13-mode reconstruction are 
obtained from an array of telescope. (One can even no¬ 
tice that for the range of angular scales considered here, 
the uncertainties obtained for an array of telescopes in 
sampling variance dominated.) 


D. Signal-to-noise ratio on r 


The computation of the signal-to-noise ratio on the 
tensor-to-scalar ratio is done by using the same Fisher 
matrix formalism as employed in the previous section, 
Eq. ([3]). For each experimental configurations and each 
value of r, we select the smallest error bars we obtained 
bin per bin. This means that for the specific case of an 
array of telescope, the estimation of the i3-mode angular 
power spectra is done by mixing different kind of sky 


apodizationtj^] 

We use the same bandpowers as in Sec. [II] and now 
add the largest angular scales, from i = 2 to £ = 20 
gathered in one single bandpower, in the summation in 
Eq. ([3]). Adding these scales is relevant for the case of 
an array of telescopes, and the case of a satellite mission. 
We will first add this bandpower at the largest scales 
for the three experimental setups, Sec. |IIID1| We will 
subsequently study its impact on the measurement of r, 
Sec. Ill D 2 (Note that we use the binned covariance 


for both the modecounting and the pure pseudospectrum 
reconstruction of Cf .) 


1. Numerical results 

Our results on the signal-to-noise ratio for r ranging 
from 0.001 to 0.2 are shown in Fig. [8] (note that for the 
specific case of a satellite mission the value r = 5 x 10~ 4 
has been added in order to fall below the 3cr limit). The 
red and black crosses correspond to a covariance matrix 
computed using the mode-counting expression for error 
bars on Cf, and the pure pseudospectrum error bars, 
respectively. The solid red line is the 3er limit. The left 
panel corresponds to a small-scale experiment covering 
~1% of the celestial sphere with a highly inhomogeneous 
noise distribution. The middle panel corresponds to an 
array of telescopes covering ~36% of the sky with a low 
level of (homogeneous) noise. Finally, the right panel is 
for a satellite mission covering ~71% of the sky with a 
low level of homogeneous noise. 

For the case of a small-scale experiment, the signal-to- 
noise ratio on r ranges from 0.06 for r = 0.001 to 4 for r = 
0.2 assuming a pure pseudo-C^ reconstruction of the 33- 
mode power spectrum (meaning that a ’’measurement” 
of r = 0.001 would be consistent with r = 0). This 
has to be compared to what would be inferred from the 
idealized mode-counting evaluation of the uncertainties, 
for which the signal-to-noise ratio varies from 0.25 for 
r = 0.001 to 6.7 for r = 0.2. For r = 0.001, 0.01 and 0.1, 
the (S/N) r derived from a mode-counting estimation of 
the uncertainties on the 33-mode is overestimated by a 
factor 4.2, 2.3 and 1.8, resp., as compared to the (S/N) r 
obtained from a pure pseudo-C^ reconstruction of the 
angular power spectrum. 

For the case of an array of telescopes and assuming 
the pure pseudosepctrum estimation of the 33-mode, the 
signal-to-noise ratio varies from 0.67 to 41 with r varying 
from 0.001 to 0.2. Values of r = 0.01 and r = 0.1 would 
be measured with a statistical significance of 5.75cr and 
28.16cr, respectively. Using instead the mode-counting 


4 We mention that mixing sky apodizations could lead to practi¬ 
cal difficulties for computing the correlations between different 
bandpowers, though this remains conceptually similar to a case 
without mixing different types of sky apodizations. 
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FIG. 8. Signal-to-noise ratio for the detection of r using pure pseudospectrum reconstruction for a potential small-scale 
experiment (/sky = 1% and inhomogeneous noise at an average level of 5.75/rK-arcminute) in the left panel, an array of 
telescope (/sky = 36% and homogeneous noise at 1/xK-arcminute) in the middle panel, and, a satellite mission (/ s k y = 71% and 
homogeneous noise at 2.2/xK-arcminute) in the right panel. Red crosses assume the mode-counting expression for the error bars 
on the reconstructed Cf and black crosses assume the pure pseudospectrum error estimated using Monte-Carlo simulations. 


estimation of the uncertainties on Cf, the (S/N) r varies 
from 3 to 54 for r ranging from 0.001 to 0.2. For the three 
selected values of r = 0.001, 0.01 and 0.1, the signal-to- 
noise ratio obtained from the mode-counting approach is 
respectively overestimated by a factor 4.5, 2.5 and 1.3, as 
compared to the realistic (S/N) r derived from the pure 
pseudo-CV estimation of Cf. 

For the case of a satellite mission, the signal-to-noise 
ratio varies from 0.66 for r = 0.0005 to 59 for r = 0.2, 
and assuming the pure pseudospectrum estimation of 
Cf. The values of r = 0.001, r = 0.01 and r = 0.1 would 
be detected with a statistical significance of 1.34, 10.84 
and 46.19, respectively. If one instead makes use of the 
mode-counting estimation of the error bars on the 13- 
mode reconstruction, the (S/N) r varies from 2 to 72 for 
values of the tensor-to-scalar ratio ranging from 0.0005 
to 0.2. For r = 0.001, 0.01 and 0.1, the mode-counting 
evaluation overestimates the signal-to-noise ratio, as 
compared to the pure pseudo-Cf reconstruction, by a 
factor 2.2, 1.59 and 1.22. 

From a qualitative viewpoint, the signal-to-noise ratios 
computed in the framework of the mode-counting expres¬ 
sion are always higher compared to the signal-to-noise ra¬ 
tios assuming the pure pseudospectrum reconstruction of 
13-mode. (This is obviously expected from the fact that 
the mode-counting approach is an idealized and under¬ 
estimated computation of the uncertainties.) We observe 
that the overestimation using the mode-counting expres¬ 
sion (as compared to the more realistic pure pseudospec¬ 
trum reconstruction of Cf ) is less marked for higher 
values of r. This behavior is common to the three ex¬ 
perimental configurations here-considered, though there 
are differences from a quantitative viewpoint. The rea¬ 
son is that for low values of r , most of the information 
comes from the largest scales, which is precisely at those 
large scales that the underestimation of the 13-mode re¬ 
construction using the mode-counting formulae is more 
marked. 


We also stress that in the case of mode-counting ap¬ 
proach, the leakages are ignored. On the contrary, 
the pure pseudospectrum approach consistently includes 
them but correct them in the analysis. This explains why 
the mode counting approach overestimate the signal-to- 
noise ratio on r. 


2. Relative importance of the reionization peak 

CMB observations covering a large fraction of the sky 
are automatically contaminated by various astrophysical 
foregrounds with complex physics involved among which 
the emission from our galaxy is the strongest. Masks are 
used to remove from the analysis the portion of sky with 
the highest foreground level, but the foreground emis¬ 
sion is present on the entire celestial sphere. Usually 
techniques - such as parametric component separation 
[5B1 used to determinate the spectral parameters or tem¬ 
plate fitting method, which deprojects the template of 
the foreground from the map EZH3H3 - are used to mini¬ 
mize the impact of the foreground. The residual level of 
foreground contaminants depends on the technique actu¬ 
ally chosen. However, the power spectrum of the galactic 
dust, polarized emission (which is the major contaminant 
of CMB measurements at frequencies above ~ 100GHz) 
behaves as £~ 2A , to be compared to £~ 2 for the CMB 33- 
mode angular power spectrum at scales above a degree 
[5]. The impact of such a galactic foreground is therefore 
expected to be more pronounced at the largest angular 
scales. 

Here, we considered the worst case scenario where the 
foreground contamination could not be removed at all 
on the largest scale, meaning that the information from 
the reionization peak is no more taken into account in 
the computation of the signal-to-noise ratio. In practice, 
we discard the first bin (2 < £ < 20) from the analysis, 
which necessarily lowers the signal-to-noise ratio on r. 
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We define this relative decrease as: 

(S/N)^ >20) - (S/N) r 
(S/N) r 


(15) 


with (S/N) r the signal-to-noise on r accounting for a/Zthe 
angular scales, and (S/N)^ >20 ' 1 the signal-to-noise ratio 
obtained by discarding the first bandpower. This relative 
decrease can alternatively be interpreted as the relative 
contribution from the first bin to the signal-to-noise on r 
since: 


(S/N)^ <20) 

(S/N) r 


(16) 


with (S/N)) £<20 ' 1 the signal-to-noise ratio on r that 
would be obtained by using the first bandpower only. 


This relative decrease of (S/N) r is shown in Fig. |j] 
The red crosses correspond to the mode-counting estima¬ 
tion of the error bars on the reconstruction of Cf while 
the black crosses correspond to the error bars from a pure 
pseudospectrum estimation of Cf. The left, middle and 
right panels respectively stand for the case of a small- 
scale experiment, an array of telescopes and a satellite 
mission. 

The case of a small-scale experiment is poorly affected 
by the removal of the first bin using the pure peudospec- 
trurn reconstruction of Cf, the relative decrease being 
systematically smaller than 0.1%. This is because in such 
a case the signal-to-noise ratio for the first bandpower, 
Cf = i/is much smaller than unity for all the 
values of r considered here. This bandpower therefore 
does not bring any significant amount of informations on 
r. This drastically differs if one uses the mode-counting 
evaluation for which <5 varies from 0.7% for r = 0.2 to 
32% for r = 0.001. This is because in this case, the 
signal-to-noise ratio in the first bandpower, Cf / 
with 2 < l < 20, becomes greater than unity though the 
sky coverage is only of ~ 1%. The fact that the relative 
decrease is more pronounced for small values of r is un¬ 
derstood as follows. For lower values of r, the recombina¬ 
tion bump at the degree scale, falls below the lensing part 
of the B -mode while the reionization bump in the first 
bandpower remains above the lensing signal. As a conse¬ 
quence, the reionization peak carries more information, 
relative to the informations carried by the recombination 
peak, for lower values of the tensor-to-scalar ratio. 

For the case of an array of telescopes, the relative de¬ 
crease ranges from 0.4% for r = 0.2 to roughly 3% for 
r = 0.001. We note that this relative decrease is now 
roughly constant from r = 0.001 to r = 0.01 and then 
decreases for higher values of the tensor-to-scalar ratio. 
This behavior of S is explained by the very same reason 
explaining why 6 decreases for higher values of r if one 
makes use of the mode-counting estimation of the uncer¬ 
tainties on the estimated Cf, and because for an array of 
telescopes, the angular power spectrum in the first band- 
power can now be measured with a signal-to-noise ratio 


greater than unity. We note that the relative decrease 
using the mode-counting behaves the same as in the case 
of a small-scale experiment (with minor quantitative dif¬ 
ferences at high values of r). 

For the case of a satellite mission, the relative decrease 
varies from 0.35% for r = 0.2 to 9% for r = 0.0005, 
with S = 8% for r = 0.001. This relative decrease 
now monotonically increases with lower values of r 
(though our results suggest that a plateau is reached 
for r < 0.001). This behavior is explained by the same 
reason explaining the decrease of <5 for higher values of r 
in the case of an array of telescopes. We also note that 
6 obtained from the mode-counting expression behaves 
the same as in the case of a small-scale experiment and 
an array of telescopes. 

As is clear from Fig. [9j the shape of S as derived 
using the mode-counting expression, is qualitatively the 
same for the three experimental configurations, though 
sky fractions and shapes of the masks drastically change. 
This is because the impact of the limited sky fraction 
is simply modelled as an overall renormalization of the 
error bars, equally applied at all angular scales (see Eq. 
(§)• Neglecting the noise contribution to the error bars 
on the .B-mode reconstruction (which is a relatively fair 
assumption here), it is easy to figure out that this over¬ 
all 1 / y// s ky does not enter in the final expression of S. 
(We note that minor differences are however expected be¬ 
cause of the different noise level and beamwidth.) At low 
values of r, the relative decrease is much less marked in 
the context of the pure pseudospectrum reconstruction of 
Cf, with respect to the mode-counting expression. This 
is because the different leakages have stronger impacts at 
large scales (in term of increase of the error bars on the 
estimated Cf), thus reducing the relative contribution 
of the first bandpower to the constraint that can be set 
on r. The impact of leakages in terms of error bars on 
Cf at large angular scales increases with smaller / sky , 
which therefore reduces the relative contribution of the 
first bandpower to the constraints on r. This is clearly 
seen in Fig. [9] where S is more important from the case of 
a small-experiment, to the case of an array of telescopes, 
to the case of a satellite mission. 


3. Performances on r detection 

As a result, at a given r and considering all the angular 
scales from 2 to 1020, the value of the signal-to-noise 
ratio is the highest in the case of a satellite mission. As 
an example, a tensor-to-scalar ratio r = 0.1 would be 
detected at a statistical significance of about 46cr. In 
the case of an array of telescope, the value of the signal- 
to-noise ratio remains high for a large range of values of 
the tensor-to-scalar ratio, showing a detection of 28cr at 
r ~ 0.1. Finally, a small scale experiment would set mild 
constraints on low values of the tensor-to-scalar ratio, 
reaching 3cr at r ~ 0.1. 
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FIG. 9. Relative decrease of the signal-to-noise ratio (in percent) if the information from the reionization peak (i.e. the first 
bin) is discarded from the analysis. The left, middle and right panels are for the small-scale experiment, an array of telescopes 
and a satellite mission, respectively. 


In the frame of the primordial B-mode detection 
prospects, the minimal value of the tensor-to-scalar ratio 
r that could be detected regarding the experimental 
setups is a relevant result. The table |TTl| summarizes in 
this perspective the aforementioned results, considering 
a measurement of r with at least a 3er statistical 
significance as a threshold. The minimal accessible 
value of r is shown with respect to the experimental 
setups and the estimation of the B-mode variance over 
all the range of multipoles (referred to as case A in 
the table). As explained above, the mode counting 
estimation of the variance overestimates the forecasts 
made on the minimal accessible r as compared to the 
realistic B-mode estimation. In the case of a potential 
satellite mission for instance, the lowest accessible r 
we could realistically expect is 2.88 greater than the 
one estimated using the mode-counting estimation. 
These results therefore highlight the inaccuracy that an 
approximative estimation of the B-mode induces on the 
performed forecasts of the detectable r values. Thus 
considering the realistic forecasts performed thanks to 
the pure estimation of the B-mode power spectrum, 
we conclude that a satellite mission would give access 
to the largest range of r, with a minimal r value of 
2.6 x 10~ 3 . A typical small scale experiment is indeed 
expected to reach only r higher than 0.1 at 3 a (note that 
r = 0.05 is detectable at 2a). Between these two cases 
lies the one of an array of telescopes, which warrants a 
detection of the tensor-to-scalar ratio if it is higher than 
5.1 x 10” 3 . As a result, each studied experiments widens 
the accessible range of the tensor-to-scalar ratio r. In 
terms of minimal detectable vlaue of r, one gains about 
a factor 20 from small-scale experiments to an array of 
telescopes, and about a factor 2 between the latter and 
a satellite mission. 

Furthermore, the table m also displays the minimum 
accessible r obtained without the information from the 
first bin (£ < 20) of the B-mode power spectrum (re¬ 
ferred to as the case B). This lack of information obvi¬ 
ously leads to a smallest range of accessible r than in the 



Small-scale 

Telescopes 

Satellite 


exp. 

array 

mission 

Mode-counting: 
Case A 

r > 0.038 

r > 0.0011 

r > 0.0009 

Case B 

r > 0.04 

r > 0.0016 

r > 0.0013 

Pure pseudo-Q : 




Case A 

r > 0.11 

r > 0.0051 

r > 0.0026 

Case B 

r > 0.11 

r > 0.0053 

r > 0.0028 


TABLE III. The minimal accessible value r with at least a 3a 
statistical significance, regarding the experimental setups and 
the estimation of the variance on the B-mode reconstruction. 
The case A means that all the bins are used. The case B 
means that the information from the reionization peak (i.e. 
the first bin from t = 2 to l = 20) has been removed. This 
is obtained by linearly interpolating the computed (S/N) r on 
our grid of values of r. 


case A for a naive estimation of the B-mode variances, 
as explained in the previous subsection. In particular, 
while the accessible r range is little affected by removing 
the first bin in the case of a small scale experiment, the 
minimal accessible r is ~ 1.44 (~ 1.04 resp.) greater for 
a satellite mission (an array of telescopes resp.) as the 
largest angular scales are relevant for these experimental 
setups. We note here that contrary to what one might 
expect, a large scale experiment would still succeed in 
detecting r of at least ICC 3 . For i between 20 and 90, 
the amplitude of the primordial signal is roughly 10% 
of the lensing signal while the total (mode-counting esti¬ 
mated) error budget varies from few percents to 10% of 
the lensing signal. Summing over the multipoles range 
thus enables a detection of r ~ ICC 3 with a 3cr statistical 
significance. 

Nonetheless, in this case B, the orders of magnitude of 
the realistic forecasts remain unchanged if the B-mode 
power spectrum is reconstructed from the pure pseudo- 
Ce approach. The values of r that could be detected at 30- 
increase by a factor of less than one percent for a small- 
scale experiment, a factor of few percents for an array 
of telescopes, and, a factor of ten percents for a satellite 
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mission. (This obvisouly reflects the values of S found 
in the previous section, Sec. HID2 ) This means that 
the pure pseudo-C/ estimation of the reionization peak 
of the B-mode mildly constraints the tensor-to-scalar ra¬ 
tio. To take full advantages of the range 2 < £ < 20 (so 
as to lower the minimal detectable value of r and to en¬ 
large the lever arm to constraint e.g. the spectral index), 
one should probably rely on more optimal techniques for 
reconstructing Cf at those largest angular scale; 


IV. CONCLUSION AND DISCUSSION 

We have investigated the detection of the tensor-to- 
scalar ratio, r, from forthcoming and potential future 
measurements of the CMB polarized anisotropies. We 
considered the B-mode as the main source of informa¬ 
tion on r and assumed the pure pseudospectrum recon¬ 
struction of its angular power spectra from the maps of 
Stokes Q and Stokes £/, previously shown to be a method 
of choice for analyzing coming data sets. We focused on 
realistic statistical uncertainties (i.e. sampling and noise 
variance) as incurred by such a numerical method, and we 
purposefully did not consider the potential gain thanks 
to delensing, nor the loss due to polarized foreground 
contamination and instrumental systematics. 

We emphasize that in this paper we consider only the 
E-to-B leakage due to a cut sky. In CMB practice there 
are numerous other potential sources of such leakages. 
For instance, they can arise from instrument limitations, 
such as beam mismatch [52| or polarimeter orientation 
uncertainty [53], or be generated by data processing, say, 
via time-domain filtering mm- Such leakages would 
also have an effect on estimated B-mode power spectrum. 
The effect will in general depend on a specific method 
used for the estimation but also on the detailed nature 
of the leakage itself, and thus would have to be studied 
cased by case. In many situations, such leakages could 
be corrected for already on the map-making stage, leav¬ 
ing therefore the cut-sky as the only fundamental source 
of the leakage to contend with on the power spectrum 
estimation level as assumed in this work. 

In contrast, we include the effects of the gravitational 
lensing, i.e., of the ’’cosmological E-to-B leakage”, in 
the total uncertainty budget, in spite of the fact that 
map-making-level, delensing procedures, which could 
correct for part of this effect have been proposed m- 
The improvements on the detection of r those methods 
could give depend on the noise level and the resolution 
of the experiment, and, on the potential use of external 
datasets (if delensing cannot be done internally). By 
including the lensing-induced B-mode, we adopt a 
more conservative viewpoint as far as forecasts on the 


5 We note that the reconstruction of CjP at large scales is also 
plagued by other sources of uncertainties such as the level of 
residual foregrounds and/or the impact of filtering of the maps. 


tensor-to-scalar ratio are concerned. 

In this framework, we first consider the case of small- 
scale (either ground-based or balloon-borne) experiments 
in an idealized way, assuming the observed sky patch 
is azimuthally symmetric. We consider four values of 
r = 0.07, 0.1, 0.15 and 0.2, and let the sky fraction to 
vary from 0.5% to 10% (with a noise level of 5.75^iK- 
arcminute at / s k y = 1%. We compare the signal-to-noise 
on r as obtained from the pure pseudospectrum recon¬ 
struction of the B-mode to the signal-to-noise ratio that 
would obtained assuming either the mode-counting esti¬ 
mation of the uncertainties on the B-mode, or a mini¬ 
mum variance, quadratic estimator. We show that the 
statistical significance on the detection of r using the 
mode-counting is overestimated by a factor ~ 1.25 as 
compared to the more realistic case of the pure pseu¬ 
dospectrum estimation. (The mode-counting also overes¬ 
timate this significance by a factor ~ 1.17 as compared to 
the minimum variance, quadratic estimator.) Similarly, 
the (S/N) r obtained from the pure pseudospectrum esti¬ 
mator is reduced by 1.5%( at r = 0.2) to 8% (at r = 0.07) 
as compared to the lossless minimum variance, quadratic 
estimator. 

For the case of small-scale experiment for which the 
reionization bump is not accessible, and in the limita¬ 
tion of azimuthally symmetric patches, the pure pseu¬ 
dospectrum approach for B-modes reconstruction is thus 
almost as accurate as the more computationally costly 
minimum variance, quadratic estimator (the former scal¬ 
ing as 0(Np{^) and the latter as 0(Np ix ) if the observed 
sky patch is not azimuthally symmetric). As shown in 
Fig. 20 of [24], non-azimuthal symmetry basically does 
not change the uncertainties on the B-mode reconstruc¬ 
tion with the pure pseudospectrum estimator (except for 
unrealistic, highly squeezed shapes). We can thus except 
this conclusion to holds for more intricate shapes of the 
observed sky. 

Our results (summarized in Tab. [I]) show that for 
a given sensitivity typical of forthcoming small-scale 
experiments, the value of the sky fraction maximizing 
the signal-to-noise ratio on r is rather insensitive to the 
method adopted to compute the uncertainties on the 
reconstructed Cf (either the mode-counting expression, 
a minimum-variance quadratic estimator or the pure 
pseudo-C^ approach). We also show that the choice of 
the bandwidth only mildly affect this optimized sky 
fraction in the case of the mode-counting approach to 
estimate the uncertainties on the B-mode reconstruction 
(see Tab. 0 - This means that using the mode-counting 
expression provide a rather reliable estimate of the 
optimized sky fraction from the viewpoint of statistical 
uncertainties though being underestimated. 

Second, we consider the detection of the tensor-to- 
scalar ratio for three selected examples, each of them 
mimicking three archetypal experimental configurations. 
Realistic sky coverage (with intricate contours) and re- 
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alistic noise distribution for the small-scale experimental 
setup are considered and the statistical uncertainties on 
the 13-mode reconstruction are derived from the mode¬ 
counting expression first (used as a benchmark), and sec¬ 
ond, from the pure pseudospectrum estimators using op¬ 
timized sky apodizations. Our results are summarized 
in the table m For each experimental setups, it shows 
the minimal values of r that could be measured with at 
least a statistical significance of 3er. One gains more than 
one order of magnitude for the minimal detectable value 
of r from the small-scale experiment to an array of tele¬ 
scopes, and another factor 2 from an array of telescopes 
to a satellite mission. This conclusion stands even if the 
largest angular scales (i < 20) cannot be used in the 
analysis. 

Let us briefly discuss the impact of those results in 
the context of single-field, slow-roll inflation. Our pur¬ 
pose here is to give a rough translation of the potential 
measurement of the tensor-to-scalar ratio with the pure 
pseudo-CV estimation of the i3-mode, into a potential dis¬ 
crimination between small fields and large fields models 
of inflation. (A more detailed study of inflationary mod¬ 
els can be found in Ref. (51], though it is restricted to 
satellite missions and assumes a different evaluation of 
the error budget for the 13-mode reconstruction.) 

The tensor-to-scalar ratio is an valuable source of infor¬ 
mation for the physics of the primordial Universe. First, 
it is a direct measure of the energy scale of inflation, 
V 1 / 4 with V the value of the inflaton potential during 
inflation: 

V 1 / 4 = 1.06 x 10 16 GeV J . (17) 

This means that a measured value of r > 0.01 corre¬ 
sponds to test a physical regime in the playground of 
Grand Unified Theories. Second, the tensor-to-scalar ra¬ 
tio is directly related to the number of e-folds, !V in f, and 
the excursion of the scalar field, A<j>, from the instant 
when cosmological fluctuations observed in the CMB are 
created during inflation, to the end of inflation [7]: 


Ai„f 



A <j> 

mpi' 


(18) 


with mpi = 1 /vo7tG the reduced Planck mass. (We note 
that N- ln { can be determined from the knowledge of the 
inflaton potential. We however let it free in order not to 
assume a too specific shape of this potential.) Single field 
inflationary models can be roughly classified between 
large fields models and small fields models, whether the 
excursion of the scalar field is transplanckian or sub- 
planckian, respectively. Though the value A(f>/mp\ = 1 
should not be considered as a sharp and univoquely de¬ 
fined frontier between small fields and large fields mod¬ 
el^ a precise measure of r then allows for discriminating 


6 The distinction between large fields and small fields models of 
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FIG. 10. Values of the excursion of the scalar field, A <f>, that 
could be observed with, at least, a 3 a significance, as func¬ 
tions of the number of e-folds during inflation. Darker blue 
to lighter blue respectively stands for small-scale experiment, 
an array of telescopes and a satellite mission. The minimal, 
detectable value of r at 3<r allowing for such a measurement 
is the one derived from the pure pseudo-C/ estimation of the 
B-mode angular power spectra. 


between this two classes of models. For TV; n f ~ 30 and 
considering zero runnings of the spectral index (see (55] 
for extensions of the Lyth bound with runnings), val¬ 
ues of r greater than ~ 0.01 would correspond to large 
fields models of inflation (see also Ref. [5S] and references 
therein for examples of small fields models evading the 
Lyth bound). 

The figure[l0 shows the ranges of A(f>/mp\ as a function 
of Ni n { accessible assuming that the tensor-to-scalar ratio 
has been measured with at least a 3cr statistical signifi¬ 
cance. Blue areas correspond to the accessible range for 
each experimental configurations (notice that the higher 
A<j>, the higher 7'). The dark blue region is for the case 
of small-scale experiments, while the somewhat lighter 
blue and light blue regions corresponds to the case of the 
array of telescopes and of the satellite mission respec¬ 
tively. The minimal detectable value of r with at least 
3cr is the one derived from a pure pseudospectrum re¬ 
construction of the angular power spectra of the i3-mode 
and using the entire set of angular scales (the case A 
of Tab. III). This shows that a measurement of r from 
the pure pseudo-C^ reconstruction of the 13-mode thanks 
to datasets coming from a small-scale experiment, is im¬ 
possible if small fields models appear to be realized in 
the early Universe. Though a detection is possible in the 
large field models, there is still a range of such models for 


inflation is usually thought in the context of effective field theory 
and the UV completion of gravity. It is therefore natural, but 
not necessary, to introduce a cut-off of the order of the Planck 
scale. 
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which the level of primordial gravity waves is still unde¬ 
tectable by a small-scale experiment. Small fields models 
are only marginally accessible from the pure pseudo-CV 
estimation of the B-mode using datasets from an array of 
telescopes, as A(f>/mp\ < 1 is accessible for N- m f smaller 
than ~ 38. A detection of r consistent with zero with a 
3cr confidence level implies an excursion of the scalar field 
(in Planck units) smaller than 0.8 to 1.8 for N - m f varying 
from 30 to 70. Finally, datasets coming from a satellite 
mission allows for a detection of primordial gravity waves 
in the small fields models with the pure pseudospectrum 
estimation of Cf, providing that the number of e-folds 
is smaller than ~ 55. On the range of e-folds considered 
here, a measurement of the tensor-to-scalar ratio consis¬ 
tent with zero then implies A(j)/mp\ < 1.2, meaning that 


a discrimination between large fields models and small 
fields models is possible for a wide range of values of 
N in{ . 
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